arXiv: 1507.04115vl [math.PR] 15 Jul 2015 


Some remarks on uniform boundary Harnack 

Principles 

M. T. Barlow* D. Karlk 
July 16, 2015 


Abstract 

We prove two versions of a boundary Harnack principle in which the constants 
do not depend on the domain. 

Keywords: Boundary Harnack principle, harmonic function, Brownian motion 
Subject Classification: 31C05, 42A61 


1 Intoduction 


The boundary Harnack principle (BHP) gives a result of the following general type. Let 
D be a domain in and £ G dD, satisfying suitable properties. Let r > 0, A 0 > 2, 
Bi = B(f,A 0 r) and B 2 = H(£,r); here £>(.,.) denote the usual Euclidean balls. Then 
there exists a constant Cd such that if u, v are positive harmonic functions on Bi fl D 
vanishing on dD fl Bi, one has 


u{x)/v{x) 

u(y)/v(y) 


< Cd for x, y G D fl B 2 . 


( 1 . 1 ) 


A BHP of this kind is called in [1] a uniform BHP, and in [13] a scale-invariant BHP. 
Here ‘uniform’ or ‘scale-invariant’ refers to the fact that the constant Cd does not depend 
on r. For Lipschitz domains D the scale invariant BHP was proved independently by 
Ancona, Dalhberg and Wu in [4] [10J ITT] . This was extended to NTA domains by Jerison 
and Kenig [12]. Bass, Burdzy and Banuelos m used probabilistic methods to obtain a 
BHP for Holder domains, but their BHP is not uniform. In [TJ a scale invariant BHP is 
proved for inner uniform domains in M d . See the papers dsns] for a further discussion 
on history of the BHP, and the various different kinds of BHP. 

In the above ‘harmonic function’ refers to functions which are harmonic with respect 
to the usual Laplacian operator in W l . (Thus these functions are harmonic with respect to 
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the infinitesimal generator of the semigroup of standard Brownian motion in M d .) Recent 
papers have studied functions which are harmonic with respect to the semigroups of more 
general diffusion type processes - see mm- 

In all these results the constant Cd depends on the domain D. For the standard 
Laplacian it is clear that such dependence is necessary, since the BHP does not hold 
for all domains D C (See however [§j where a result of BHP type with constants 
independent of the domain is proved for harmonic functions with respect to fractional 
Laplacians.) 

This paper originates in the work of Masson [lp] , where a certain kind of BHP with 
constant Cd not depending on D was needed - see m Proposition 3.5]. Masson’s work 
was in the context of discrete potential theory for Z 2 . Let S x = (S%, k > 0) be the simple 
random walk on Z 2 , started at x, and write S = S°. Write Z 2 , = {(op, x 2 ) G Z 2 : x\ < 0}, 
and let Q(x,n) = {y G Z 2 : \x — y\ < n }. Let N > 1, K C Q(0 : N) D Z 2 , and 
D = Q(0, N) — K. (The case of interest is when 0 G K.) Let r + = minjfc > 1 : Sk ^ D}, 
and F = {S T + G Z 2 — K}, so that F is the event that S leaves Q(0, N) before hitting K. 
Let W = {(xi,x 2 ) : 0 < |x 2 | < aq}, so that IF is a cone with vertex (0,0) and angle 7 t/ 4. 
Massons’s theorem is that there exists po > 0, independent of N and K, such that 

ns r+ G W\F) > p 0 . (1.2) 

This result extends to give also P(5'* + G W\F) > p\ for x = (xi,x 2 ) G Q(0, N/16) with 
X\ > 0. The fact that Pi do not depend on the structure of K is essential in the context 
of |T5], since K is a random path (actually a loop erased random walk), and the estimate 
(1 1.2 j) was needed for all possible K. 

Although the connection with BHP is not made in [15], this result is clearly of BHP 
type. For x G Z 2 let r = min{fc > 0 : Sk ^ D}, and define the functions 

n(x) = P(S'* G K c ), u(x) = P (S x G K c D W). 

These are (discrete) harmonic in D , and P (S x G W\F) = u(x)/v(x). Since u < v it is 
immediate that 

<Pi\ for x, |/ G Q(0, N/16) fl (Z 2 — T?_). (1.3) 

n{y)/v{y) 

Thus we have a BHP in which the constant Cd = Pi 1 does not depend on A"; the price 
is that the inequality only holds for x close to 0 with X\ > 0. 

In this paper we present two further results on this kind. In Section [2] we extend 
Masson’s result to Euclidean space, and prove it for d > 2. (The result is trivial for 
d — 1). In Section [3] we prove a BHP for d = 2 where the constant Cd does not depend 
locally on K, but does require some regularity of K away from 0. A counterexample 
shows this result does not extend to d > 3. 

Throughout this paper we write X = (X t , t G [0, oo), P x , x G M d ) for Brownian motion 
in P x is the law of X started at x. For a set A C we define 

Tfi = inf(F > 0 : X t G A }, ta = = inf{t > 0 : X t 0 A}. 
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We write B(x,r) = {y : \x — y\ < r} for Euclidean balls. When we use notation such as 
C = C(a, d ) this will mean that the (positive) constant C depends only on the parameters 
a and d. 


2 Uniform BHP for a harmonic functions associated 
with cones 

In this section we will prove the Boundary Harnack Principle for two fundamental har¬ 
monic functions. 

Let 

H_ = {x = (xi,..., Xd) G R 0 ' : X\ < 0}, B~ = 5(0,1) D EL, 

and define H + , B + analogously. We write B = 5(0,1). Let K C EL be closed and 
connected. Set 

D = 5(0,1) - K. 

Let TT\ : R^ —^ R d be projection onto the ay-axis, so 7r((xi, x 2 , ■ ■ ■, Xd)) = (xi, 0,..., 0). 
Let W a be the cone 


W a = {z G R d : \z — 7Ti(z)| < z\ tan(a)}. 


Set 

W a {r) = 5(0, r) D W a . 

Write t — Td for the exit time of X from D, and define the functions 

v(x) = F x (X t edD n K c ), u(x) = P- r (X T edD n w a ). (2.1) 

Thus u, v are bounded, positive harmonic functions which vanish on dK, and have bound¬ 
ary values 1 on dD D K c and dD D W a respectively. It is clear that u < v on D. Both 
u and v are bounded, positive and harmonic inside the domain D. Hence they satisfy 
the usual Harnack Inequality, [6], Theorem II. 1.19], in balls which are far enough from 
the boundary of D. The main result of this section is that these two functions satisfy 
a BHP with constant which depends only on d. (Note that since the geometry of the 
boundary of K is not specified, classical results on the Boundary Harnack Principle such 
as P, Theorem III. 1.2], do not apply.) 

The main result of this section is the following Theorem. 

Theorem 2.1. Let a G (0,7r/2]. There is a constant C = C(a,d ) > 0 depending only on 
a and d, and independent of K, such that for any x,y G 5(0,1/2) D EI + , 

LLhd < c . ( 2 .2) 

v{x)/v{y) 

Remark 2.2. (1) The result does not require 0 G K, though this is the most delicate 
case. 
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(2) Let x = (1/4,0). Then the usual Harnack inequality gives that u{x) > C'(a,d), and 
so, since v < 1, (12. 2 p implies that 

C'(a, d) < u(y)/v(y) < 1 for y G 5(0,1/2) D H + . (2.3) 

On the other hand, given (12.3(1 the inequality (12. 2 p follows immediately. 

(3) Suppose that (12.211 (and therefore (12.311 ) holds for some /3 G (0,7t/2), and let a G 
(/3,7r/2]. Writing u Q , up for the functions given by (12.ip we have up < u a < v < 1. So 
([23]) for up implies (12.3p . and therefore Theorem 12.11 for u a . ft is therefore sufficient to 
prove Theorem 12.11 for j3 G (0,0^), where the angle depends on d. 

Let Kd = [4(d + 2)] -1//2 , and define ad = arcsin^) for d > 2 . Set 

A ad = w ad ndB(o,i). 

Note that the proportion of this surface A ad to the surface of the unit ball decreases as 
dimension grows. Similarly, let f3 G (0, ad) and define Kp = sin(/3) and 

Ap = WpAdB{ 0,1). 

First we observe that by Lemma 12.31 below it is enough to prove this statement for 
7Ti(x) and 7Ti (y). 

Lemma 2.3. If x G W ad (l/2) then 7ii(x) G W ad ( 1/2), and there is c = c(ad) such that 
c^ 1 f (tti(x)) < f(x) < c/(7Ti(x)) for any non-negative bounded harmonic function f in D. 

Proof. By definition of the cone, it is clear that -Ki(x) G W ad { 1/2). If d > 3 then x G 
B(tti(x), Xi tan(ct d )) C B(tti(x),x i) C B + , and applying Harnack’s inequality, the result 
follows. If d — 2 then it is enough to consider the Harnack chain with two balls and to 
apply the Harnack inequality twice. □ 

We now consider the exit distribution of Brownian motion from the unit ball B , 
without considering the set K. Set 

hp(x) = ¥ x (X TB eWp). 

This is a harmonic function on B, and potential analysis provides a density function which 
is the d- dimensional Poisson kernel 

Pd{x,y) = w~ d \]^ J^L; 

here Wd-i is the surface area of the unit sphere with respect to the surface measure cTd{dy ) 
on the d-dimensional unit sphere. Hence hp admits an integral representation given by 

hp(x) = / P d {x, y)a d {dy). 

JWpndB 

(See j6j Theorem II.1.17]). 

Our argument is based on two main steps: comparison of the values of hp inside the 
ball and the connection of hp with two functions u and v. For this purpose, we compare 
first the values of hp on the left-half of the ball, B~, with the values on the positive axes 

{(xi,0,..., 0) : 0 < xi < 1}. 
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Lemma 2.4. For any x = (aq, 0,..., 0) with X\ G (0,1) 

MO) < hp(x). 

Proof. Define the distance function 


d(z,E) = M{\z-y\ :y G E}. 

Set A 1 = dBDdW ■ thus Ai is a d— 2 dimensional sphere. Let rq = d(x, Af)\ by symmetry 
x is the same distance from all points in A 1 . Set 


B x = B(x,r i). 


This is the ball with center at x whose surface crosses the surface of the unit ball 
through A l . Then A p = dB n B x , and writing o' d {-) is the surface measure on dB x , 

-B) ai {A t ) 

r,j(dlP) ~ o d (dB )' 

Note that a Brownian motion X started at x and leaving the ball B x through dB x — B 
must leave the unit ball through Ay. So 


hy{x) = W{X TB e Ay) > W{X TBX G dB x - 
a' d (dB x -B) o d (Ay) 
a ' d {dB x ) - a d (dB ) 1 


B) 


G Ay) — hy{ 0). 


□ 


The final piece of argument is based on the comparison of the values of hy. For this 
purpose, we need the following two technical lemmas. 


Lemma 2.5. Let d > 2. There is a dimensional constant n d 
x G [0, K d \ 


Kd - x 

K d + X 


1 + x 2 — 2xn d 

1 + X 2 + 2xK d 



1 


so that for any 


Proof. Denote the functions on the left hand side and the right hand side of the inequality 
by f(x) and g(x), respectively. It is clear that /(0) = g(0). Hence it is enough to show 
that for every x G [0, K d \ the derivatives satisfy f'{x) < g'(x). Here 


/'(*) 

g\x) 


(«d + x ) 2 

2 n d (d + 2)(1 — x 2 )(l + x 2 — 2xK d ) d l 2 
(1 + X 2 + 2xK d ) 2+d / 2 


Now using 0 < x < n d , it is not difficult to see that 


-g'{x) < 2n d (d + 2) 


2 n d 2iK d 

{2n d ) 2 ~ (K d + x) 2 


~f\x). 


This inequality together with the Mean Value Theorem leads to our result. 


□ 
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Lemma 2.6. Suppose K d is as defined before. Then for any x, z G [0, K d ] 


f(z,z) 


1 — x 2 1 — x 2 

(1 + x 2 — 2 xz) d / 2 (1 + x 2 + 2 xz) d / 2 


< 2 . 


Proof. The result follows immediately for x = 0. First we fix x G (0, K d \ and consider the 
change in the direction of z. It is easy to show that jpf(x,z) > 0. Hence the maximum 
occurs at z = K d , that is, f(x,z ) < f(x,K d ) for any x, z G [0, K d ]. Next we differentiate 
f(x, Kd) with respect to x. 

d T d(l - x 2 )(K d - x) d(l - x 2 )(/c d + x) 

Q x {X ’ M ~ [(1 + x 2 - 2xK d ) l+d / 2 (1 + a; 2 + 2xK d ) 1+d / 2 _ 

2x 2x 

(1 + x 2 — 2xK d ) d l 2 (1 + x 2 + 2xK d ) d / 2 _ 

The difference inside the first parenthesis is negative by Lemma [2751 Therefore the deriva¬ 
tive is negative and the function /(•, K d ) reaches its maximum at x = 0 where /( 0, K d ) < 2. 
So f(x, z ) < f(x, K d ) < /(0, K d ) < 2 for any x,z e [0, K d ]. □ 


The previous lemma helps us to prove the following statement. 


Lemma 2.7. Let d > 2. Then for any x e B , hg{x) < hp( 0). 

Proof. Since hp is harmonic in B and continuous on B~, it reaches its maximum on 
dB~ = (dB HdB~) U ({(xi, x 2 ,..., x d ) : x x = 0} C\dB~). Since hp is zero on (dB DdB~), 
it is enough to find the maximum of hp on the set {(xi, x 2 ,..., x d ) : x x = 0} D dB~ and 
to show that this maximum value is bounded above by hp(0). 

We can reduce the set of interest to an even further subset. The rotational invariance 
of Brownian motion together with the symmetry of the domain give that hp(x) = hp(y) 
for any x — (0, x 2 ,..., x d ) and y — (0 ,y 2 ,..., y d ) with \x\ = \y\. Hence we only need to 
consider the points of the form x = (0, x 2 , 0,..., 0) with x 2 G [0,1). 


First note that for any z — (z x ,..., z d ) G Ap, \z 2 \ < Kp < K d . So if Kp < x 2 < 1 then 

- o d (dz) < 0. 


(2.4) 


dhp 1 

f -2x 2 | 

x — z\ 

2 - d(x 2 - z 2 )( 1 - x\) 

dx 2 C Od-! J 

Aa d 

x — z\ 

\d+2 


Hence hp is a decreasing function in x 2 whenever x 2 G ( Kp, 1) and 

M(0, X 2 , 0,..., 0)) < ^((0, Kp, 0, ..., 0)) 


for Kp < x 2 < 1 . 


Assume that 0 < x 2 < Kp < K d . Split the surface Ap into 2 n parts as follows: Take 
n+1 non-negative numbers {C(f)}^ =0 such that 0 = £(0) < £(1) < £(2) < • • • < ((n) = Kp 
and define strips {5' J }" =I in a way that 

Si = {z = (z x ,..., z d ) G Ap : z 2 G [C(i - 1), C(*))} 
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and the measure of each strip, a d (Si), equals each other. Similarly, define n negative 
numbers {£(i)}^l_ n by £(—*) — — C{i) f° r * = 1, Also define the strips, {S'_j}”T 0 1 , on 
the lower-half of Ap the same way as above 

S-i = {z = (z u . . ., z d ) e Ap : z 2 e [C(—* - 1), C(—*))}■ 

Then 

Ap = IJ Si and a d (Si) = a d (S j ) = i, j e {-n + 1,..., n}. 

—n+l<i<n 

Note that if z — (zi ,..., z d ) £ Si then 

\x — zj 2 = 1 + x 2 — 2x 2 z 2 > 1 + x 2 — 2x 2 ( (i). 


Using this partition 

hp(x) = / P d (x,z)a(dz) = ^ / p d(x, z)a(dz) 


Z—— n~\~l' 


n —1 


P d (x,z)a(dz) + / P d (x, z)a(dz) + ^ 


I So 


i= l 


P d {x : z)cr(dz) + / P d (x, z)a(dz) 
Js-i 


< 


i-x\ 


u d -i L(1 + x\ - 2x 2 K d ) d / 2 

n —1 r i 2 

1 — xi 


cr d {S n ) + — — — ^rjjz^d(So) 


E 


uj d _i L(! + x \ - 2z 2 CW) d/2 


(l + x 2 ) d / 2 
&d{Si) + 


l-x\ 


(1 + - 2x 2 ( (-i)) d P 


Vd(S_i 


< 


1 ~x\ 


+ 


i-x\ 


(?d(Ap) 

uj d -i -2n L(l + x 2 - 2x 2 K d ) d / 2 ' (1 + xl) d / 2 
a d (Ap) "- 1 r 


uj d -i ■ 2n 


E 

i=l 


l-x 2 


+ 


l-x\ 


(1 + x\ — 2x 2 ( (i)) d P (1 + x\ — 2x 2 ((—i)) d / 2 


By the restriction 0 < x 2 < Kp < n d , the first term is bounded by c d /2n where c d = 
2a d (A ad )/u d -i(l — 2n 2 d ) d l 2 . For the term inside the second bracket, Lemma 12761 provides 
an upper bound and hence 


x% 


+ 


Xn 


(1 + xl - 2x 2 ( (i)) d / 2 {l + x 2 2 + 2x 2 ( (i)) d / 2 


< 2 . 


Hence for any n e Z + we obtain 




Qd{Ap) n- 1 
Ud-i n 
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(2.5) 


Finally, if we take the limit as n —> oo then 

hp(x) < = hp( 0). 

Ud-i 

for any x = (0, ay, 0..., 0) with ay G [0, up). 

By (12.4p and (12.51) . we conclude that hp(x) < hp( 0) for any x G B~. □ 

Proof of Theorem \2.1\ Let r G (0,1], 0 < /3 < ad, and write h for hp. Let u be the 
function defined in (12.ip with angle (3. 

First, let x = (ay, 0,..., 0), with 0 < ay < 1/2. By the Markov property, 

W[h(X TB )\T K < t b ] = E XT *[h(X TB )\ = h(X TK ) < sup h(y). 

y£K 

By Lemma [2.41 and Lemma [2.71 we have h(y) < h( 0) < h(x) for any y G K. Thus 

E x [h(X TB )\T K <r B \<h(x). 


So 


h(x) = E x [h(X TB )] 

= W[h{X TB )\T K < t b \ • P*[T* < t b ] + E x [h{X TB )- t b < T k \ 
< h(x)( 1 — v(x)) + u{x). 


Hence 


Cfl = h( 0) < h(x) < 4^ = V[X T , e A„\t b < T k ], 


(2.6) 


This proves (12. 2 p for x on the ay-axis, and by Lemma 12.31 it then follows for x G Wp. 

Now let x G i?(0,1/4) D H + , and set x' = x — tti(x). Let W' — x' + Wp, A' = 
W' ndB(x', 1/4), and write r' = t B (x',i/4)- Then applying (12.6p to the ball B(x', 1/4) — K, 
we obtain 

cp < P' r [X r , G Ay < T k \ = f f' ,T < Tk \ (2.7) 

p ~ 1 J F x [t'<T k \ k j 

Since r' < r = t Bi this implies 

W X [X T , G A'; r' < T k ] > cpv(x). 

Now by the standard Harnack inequality, 

¥ y [X r <E Ap,r < T k ] > ci for y G A'. 


u{x) = W[X T eAp,r<T K ] 

> P"[^r G Ap, X T i G A!, r' <T k ,t< T k ] 


= E" 


<t k ,x t ,&a’)^ X tI [X T G Ap,r < T k ] 


> ciP x [t / < T k ,X t i g A!] > cicpv(x). 


Then 


( 2 . 8 ) 

















Since we have u(x) < v(x) everywhere, it follows that 

C(8,d) < < i for x e e+ n b( o, 1/2). (2.9) 

v{x) 

This proves (j2.3[) for /3, and Theorem 12.II then follows from Remark 12.21 □ 


3 General Case: BHP for positive harmonic func¬ 
tions in d = 2 


In the previous section, we proved a uniform Boundary Harnack Principle for a particular 
pair of harmonic functions. In this section, we give a second result of this type, where the 
constant in the BHP does not depend on the structure of dD near the boundary point 0. 
The result is proved for d = 2 only, and a simple counterexample which is given at the 
end of the section shows it does not hold for higher dimensions. 

We define the quasihyperbolic metric kn(x,y) for a bounded domain h C R 2 by 

kn(x,y) = inf / 

7 d(z , <9fi) 


where the infimum is taken over all rectifiable curves 7 connecting x to y in and d(x, dfl) 
denotes the distance of x to the boundary of the domain fi. It should be noted that 
the metric kn is closely related with Harnack inequality for positive harmonic functions. 
To see this, we say {5(xj, 7 y )}™ =1 forms a Harnack chain connecting x G 12 to y G 
if B(xi,2ri) C H, x G B(x 1 ,r 1 ), y G B(x n ,r n ) and B(x h ri) n B(x i+ i, r i+1 ) ^ 0 for 
% — 1,..., n— 1. The shortest length of a Harnack chain connecting x to y is comparable to 
kn(x, y) +1. Thus we can rewrite the Harnack inequality for a positive bounded harmonic 
function h on the domain as 


exp (-c(ka(x,y) + 1)) < 


h{x) 

Ky) 


< exp (~c(k n (x,y) + 1)), 


x,y E fl. 


Definition 3.1. We say that satisfies the quasihyperbolic boundary condition if there 
exists Xq G and constants C\ and C 2 such that 


kn(x, xo) < Ci log(d(x 0 , dfl)/d(x, 90)) + C 2 , x G fl (3.1) 

We will make us of the following two lemmas concerning quasihyperbolic metrics. 

Lemma 3.2 (Smith-Stegenga [TSj). Let xq G and C 5(0,2) satisfy the quasihyper¬ 
bolic boundary condition with constants C \, GV Then there exist constants r > 0 and c\, 
depending only on d, C\ and C 2 , such that 


exp(r kn(x, xo))dx < c\. 


(3.2) 
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Ill [IS] it is only stated that the integral in (13. 2 p is finite; however the proof there 
gives that it is bounded by a constant depending only on the constants Ci,C 2 in (13.111 , 
the dimension d, and the diameter of 

Lemma 3.3. Let u be a non-negative subharmonic junction on a bounded domain C M 2 . 
Suppose there is 5 > 0 such that 

Is := / (log + u(x)) 1+s dx < oo. 

Jn 


Then 


u(x) < exp(2 + ciy s d(x, dLl) 2 ^ 5 ), 


where c > 0 is a constant depending only on 5. 


While the main argument of this lemma is due to Domar m the proof of this lemma 
as stated here can be found in [3] Lemma 3.1]. 

For the simplicity of notation, let us denote the unit disk and the disk with the center 
at (1,0) and radius 1 by B 0 and Bi, respectively. Consider a compact and connected set 
K C B 0 — Bi. For r G (0,1) and e G (0,1) we define two hitting angles 


6i(r,e,K) = inf{# G (0,27t) : B{re ie ,re) OK 0}, 
9 2 (r, e, K) = sup{0 G (—27 r, 0) : B(re ie , re) D K ^ 0}, 


and a tube around the arc from 0 2 to 0\ 

F(r, e, K) = ]J B(re ie ,re) 

e 2 (r,e,K)<0<0 1 (r,e,K) 

assuming there exists at least one 6 G (0, 2tt) such that B(re ie ,re) D K ^ 0. 

Definition 3.4. We call a set K c B 0 — Bi as (r, e)-good if there are two points x\ G 
(dK)C]B(re iei ,re) and x 2 G (dK) C\B (re id2 , re) such that B{x\,re)—K and B(x 2 ,re) — K 
satisfy the quasihyperbolic boundary condition. These two points x\ and x 2 will be called 
roots of K. 

In this section, we consider on an (r, e)-good set K with roots x± and x 2 , and with 
constants C\ and C 2 used in (13.If) . Unless otherwise stated all constants c: t will depend 
only on C\,C 2 ,d and e. 

The main result of this section is the following theorem. 

Theorem 3.5. Suppose K C B 0 — B i is a compact and connected, ( r,e)-good set for 
some | < r < | and e < |. Suppose that u and v are positive, bounded and harmonic 
on D = B 0 — K, continuous on D and vanishing on (dD) D Bq. Then there is a positive 
constant c = c(e, d, C±, C 2 ) such that 

u{x)/v{x) < forx,y G 5(0,1/32) DD. 
u{y)/v{y) 
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We begin by proving a rooted local Carleson estimate. 

Lemma 3.6 (Carleson estimate). Suppose u and K are as in Theorem, \S. 51 and x\, X 2 
are two roots of K corresponding to 9\ and 62 , respectively. Let wo = (r, 0) G B 0 . Then 

u(x) < c £tCliC 2 u(w, 0 ), x E B{xi, re/2) — K, i = 1, 2. 

Proof. We will prove this for x \; the same argument also applies to x 2 - First denote 
B(re ldl ,re) U B(xi,re) by B. By our assumption, u is a positive and bounded harmonic 
function in B — K vanishing on (dK) D B. Extend u to B by setting it equal to zero on 
K D B and take its upper-semicontinuous regularization. Without loss of generality we 
denote this extension by u as well. This function u is non-negative sub-harmonic in B 
and harmonic in B — K. Then 



Since u is harmonic in B — K, we obtain using the Harnack inequality 

u{l%) ~ C ex P ( ck B- K (x,re id1 )). (3.3) 

Now we need the following elementary inequality: given positive constants a > 0 and 

log ” w - (iyy) =(i) xB ' x - h 

To see this, it is enough to observe that derivative of log(x) is bounded above by the 
derivative of x^^ a /{(5/a) for x > 1. Making suitable choices for a and /?, and using this 
inequality with (13. 3 j) and Lemma 13.21 gives 


h < c 


exp (t k s _ K (x,re l<>1 )) < ci, 


' B-K 


Now if we use Lemma [3731 and restrict x to B(xi,re/2) 
and hence 


u(x) 
u(re ldl ) 


< c 


e,ci,C2- 


K then d(x,dB) > re/2 > e/8 


Finally we can End a Harnack chain in F(r,e,K ) connecting two points re lSl and Wq to 
show that u(re ldl ) < c e}Cl}C 2 u(wo). The result now follows. □ 


Proof of Theorem, \3.FA First we define the following three paths which enclose the region 
we will work on. 


Ei{r) \= {tx 1 + (1 - t) re 191 : t e (0,1)}, 
E 2 (r) := { tx 2 + (1 - t)re 102 : t G (0,1)}, 
E 3 (r) := {re 19 : 0 2 (r,e,K) < 9 < fo(r,e, W)}. 
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Let 0(r, e, K) be the domain enclosed by the curves Ei(r), E 2 (r), E 3 {r) and the and the 
set K. Define E±{r) = 90(r, e, K) — (5i(r) U E 2 (r) U E 3 (r)). Since u is harmonic in D , 
for any x G [5(0,1/32) D D\ C 0(r, e, K) 

u(x) = [ u{y)¥ x \X Tn(re K) G dy) 

JdO(r,e,K) 

= / «(y)r(X ra(r . K| 6 dy). 

J E\(r)\jE2{r)\jEz(r) 

By the Carleson estimate Lemma [3761 u is bounded above by c £jCljC 2 u(w 0 ) on 5 1 (r)U5 2 (r). 
By using a Harnack chain as in the proof of Lemma 13.61 and the regular Harnack inequality, 
we also obtain that u(y) < c £tC 1 }C 2 u(w 0 ) for y G 5 3 (r). Hence 

u(x) < c £iCliC2 u(w 0 ) ¥ x (X Tu(r ^ K) G 5i(r) U E 2 (r) U E 3 (r)). (3.4) 

Second, we define a tube F around 5i(r) U E 2 [r) U 5 3 (r). 

F= (J B(y, re/ 4). 

y&E\ (r)VjE2(r)\jEz(r) 

Let ip be smooth function with compact support in F such that ip — 1 on E\ (r) U E 2 (r) U 
E 3 (r). We can choose ip so that \Aip\ < c d . Then for x G 5(0,1/32) D D 


[ *P(y) px (^r n( ., e , x) e dy) = ip(x) + [ a ip(y)G a ^ e>K )(x, y) dy. 

J dfl(r,e,K) J Q(r,e,K)r\supp(tjj) 

(see [2, Lemma 1] for details). This equation yields to 

px (^- n(reK) e5 1 (r)U5 2 (r)U5 3 (r)) < [ \AiP(y)\G D (x,y)dy. (3.5) 

J fl(r,e,K)nF 

Gd{x,-) is a positive harmonic function in the domain D — 5(0,1/32) containing 
0(r, e, K) D F. Moreover, 

<9(0(r, e, K) n F) = [9(0(r, 6, K) n F) D dD] U [9(0(r, e, K) n F) n D] 

= [9(0(r, e, K)HF)n dD) U [90(r, e, K) n D] U [dF n 0(r, 6, K)\. 

Lemma [3761 and the Harnack inequality apply to the function Gd(x, •) on 90(r, e, K ) fl D 
so that 

G d (x, y) < c £ , CuC2 G d (x, w 0 ), ye 90 (r, 6, K) ft D. (3.6) 

as in the first part of the proof. Next, 

dF n 0(r, e, K) = [dF D 0(r, e, K) n 5(xi, re/2)] 

U [dF D 0(r, e, K ) fl B(x 2 , re/2)] 

U [dF fl 0(r, e, K) — (B(x i, re/2) U B(x 2 , re/2))]. 
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We note that any point on the third part is at least at a positive distance away from the 
boundary of K and this distance is bounded below by a constant depending on e. On 
the first two parts, we can use Lemma 13.61 and on the third part we can use the Harnack 
inequality on a Harnack chain to prove that 

G d (x, y) < c e>CUC2 G D (x, w 0 ), yedFD 0(r, e, K). (3.7) 

So by (13.61) and (j3.7|) we obtain 

G d (x, y) < c e , cljC2 G D (x, w 0 ), ye d(F n fl(r, e, K)), 

and by maximum principle this inequality also holds inside the domain F D Q(r,e,K). 
Then Inequality (13.51) together with the last inequality show 


pX (^r n(7 . e , x) 6 E^r) U E 2 (r) U E 3 (r)) < c £ jC1jC2 c'G D (x,w 0 ), (3.8) 

where d = f n(r ^ K)n p \^(v)\ d V- 

For the final part of the proof, denote the point (r/2, 0) by wq and consider the circle 
dB(u>o,r /8). By our assumption on the harmonic function v and the Harnack inequality 


Moreover 

and so 


v(z) > cv(w 0 ), 

G d (z, wq) < c, 
v(z) 


G d (z, w 0 ) < c 


z G dB(w 0: r/8). 
z G dB(w 0 ,r/8) 
z G d(D — B(w 0: r/ 8)) 


v{w 0 y 

since v is positive. By the maximum principle, the last inequality holds inside the domain 
D — B(wo,r/8) which includes B( 0,1/32) D D. Using this inequality, 
obtain 

v(x) 


and 


we 


u(x) < c e 


/ \ v(x) 

u(Wq) < c, 


u(w 0 ) 


v(w 0 ) ~ £ ’ C1 ’ C2 v u/ V («; 0 ) 

where we applied the Harnack inequality to compare v(w 0 ) to v(w 0 ). 

Finally, if we switch the roles of u and v and the roles of x and y we also obtain 


which leads to the result 


v(y) < v i w o) 

u(y) ~ Ce ’ Cl ' C2 u(w 0 ) 


u{x)/v{x) 
u(y)/v(y) ~ 


□ 
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We conclude this paper with an example which shows that one cannot expect a uniform 
BHP in higher dimensions, even if the set K is contained in EL and we restrict the 
inequality to x,y in {x : 7Ti(x) > 0}. Let B = B( 0,1), K 0 = B D Ho, where H 0 = {x : 
n\{x) = 0}. Let e, 5 be small and positive. Set 

K = K 0 - B( 0, 5), D = B( 0,1) - K. 

Let y be on the x\ axis with TTi(y) = 1/4. Let w_ and u + be the harmonic functions in 
D with boundary condition 1 on dB D H_ and dB D H + respectively, and zero boundary 
conditions elsewhere. Set v = + u+. So if r = To then we can write 

u_(x) = W(X T G dD D EL, r < T K ), v(x) = P x (r < T K ). 


By symmetry we have 


On the other hand if B' 


Mo) 

w(0) 

B( 0, 5) then 


1 / 2 . 


P y (T B , < t d ) < P y {T B , < t b ) < c8 d ~ 2 . 


So we have 

v(y) x 1 , u-(y) < c8 d ~ 2 . 


Thus 


M-(0)/n(0) 2 _ d 

u-{y)/v{y) - 


(3.9) 


By continuity this inequality will also hold if 0 is replaced by a point x close to 0 with 

Ti(x) > 0 . 
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